APS/123-QED 



Spatial structure of Mn-Mn acceptor pairs in GaAs 
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The local density of states of Mn-Mn pairs in GaAs is mapped with cross-sectional scanning 
tunneling microscopy and compared with theoretical calculations based on envelope-function and 
tight-binding models. These measurements and calculations show that the crosslike shape of the 
Mn-acceptor wave-function in GaAs persists even at very short Mn-Mn spatial separations. The 
resilience of the Mn-acceptor wave-function to high doping levels suggests that ferromagnetism 
in GaMnAs is strongly influenced by impurity-band formation. The envelope-function and tight- 
binding models predict similarly anisotropic overlaps of the Mn wave- functions for Mn-Mn pairs. 
This anisotropy implies differing Curie temperatures for Mn £-doped layers grown on differently 
oriented substrates. 

PACS numbers: 71.55.Eq, 73.20.-r, 75.50.Pp, 75.30.Hx 



The properties of dilute ferromagnetic semiconductors, 
such as Gai_ x Mn x As, depend sensitively on the nature of 
the spin-polarized holes introduced into the host by the 
magnetic dopants Q, 0. Considerable controversy per- 
sists about the nature of isolated magnetic dopants in 
many semiconducting hosts (e.g. Mn dopants in GaN). 
Measurements of the local density of states (LDOS) near 
an individual Mn substituted for a Ga atom in GaAs 
(MnQa) by cross-sectional scanning tunneling microscopy 
(X-STM) 3] have resolved this question for Gai_ x Mn x As: 
there is a hole state bound to the Mn dopant, yielding a 
Mn 2+ 3<i 5 + hole complex jj, y, [y, [j| that produces an ex- 
tended, highly anisotropic LDOS. The anisotropic shape 
of the bound hole state at distances > 1 nm, originat- 
ing from the cubic symmetry of GaAs suggests highly 
anisotropic Mn-Mn interactions}^. [ToL ll l| . 

The ferromagnetic properties of Gai_ x Mn x As, how- 
ever, depend on whether this shape persists for con- 
centrations x of Mn impurities for which Gai-^Mn^As 
is ferromagnetic (x > 0.01) [l2j. Popular models of 
Gai-^Mn^As assume that holes residing in a bulk GaAs- 
like valence band, and thus evenly distributed through- 
out the material, mediate the ferromagnetic interaction 
among Mn spins Il2l. However, angle resolved photoemis- 
sion spectroscopy [l3| observes an impurity band near Ep 
and infrared absorption measurements reveal a strong 
resonance near the energy of the Mn acceptor level as 
well as deeper in the band-gap of GaAsfTS. HEj . Fur- 
thermore, recent Raman scattering experiments have 
shown that a Mn 2+ 3d 4 configuration partially occurs for 
x > 0.02 16]. Recent theories suggest significant modifi- 
cations in the ferromagnetic properties of Gai-^Mn^As 
if the holes reside in a strongly disordered impurity 
band|Tll 0|i|i|23. If the Mn density is near the 
metal-insulator transition, and individual Mn dopants 
states are weakly hybridizing with each other, then the 
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FIG. 1: (Color online) Room-temperature X-STM constant- 
current image of a section of a Mn 8- doped layer with inten- 
tional Mn concentrations of 3 x 10 13 cm -2 : (a) 51 x 29 nm 2 ; 
(b) 12 x 10 nm 2 ; (c) 13 x 10 nm . Images were acquired at a 
sample bias of U s = +1.5 V on a cleaved (110) surface. 



inhomogeneous hole density of the impurity band mea- 
sured by X-STM near an individual Mn dopant should 
closely resemble that of an isolated neutral Mn, only 
weakly perturbed by Mn-Mn interactions. 

Here we present experimental evidence that the shape 
of a Mn acceptor in a Mn-Mn pair remains anisotropic 
and retains the crosslike shape of a single Mn even 
when the dopants are separated by only 0.8 nm, which 
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FIG. 2: (Color online) Topography of Mn-Mn pairs: (a) 12 x 
12 nm 2 , X-STM image of two ionized Mn separated by 1.4 nm 
(Us = —0.6 V); (b) Schematic of the (110) surface showing 
the location of subsurface Mn; (c) same area as (b), imaged 
with U s = +1.1 V so the Mn are neutral; (d) 6 x 6 nm 2 , X- 
STM image of two Mn separated by 0.8 nm (U s = +1.55 V)\ 
(e,f) Calculation with TBM of (c,d), with parallel spins; 
(g,h) Difference between Mn-Mn pairs with parallel spins and 
non-interacting single Mn's; (i,j) Same as (g,h) but antipar- 
allel spins. 



is the typical distance in Gao.96Mno.04As. The over- 
lapping Mn wave-functions at such short Mn-Mn sep- 
arations are exceptionally hard to disentangle in bulk 

GaMnAsHJill^. In the Mn ^ do P ed la ^ ers of Fi g- □ 
however, isolated pairs and dense groups of Mn accep- 
tors at these distances can be analyzed separately from 
surrounding dopants. 

The measurements were performed on several samples 
using chemically etched tungsten tips. The J-doped lay- 
ers were grown at 370 °C by molecular beam epitaxy 
on a 100 nm GaAs buffer on a Zn-doped (001) GaAs 
substrate. The high growth temperature was chosen 
to suppress the appearance of structural defects such 
as As antisites, which would complicate Mn identifica- 
tion. The higher growth temperature also led to in- 
creased segregation, which broadened the £-doped layers 
of Fig. H Despite the high growth temperature a low 
density (< 10 17 cm -3 ) of As antisites was observed and 
clearly identified as charged n-type donors (not shown in 
the figure). The J-doped layers themselves clearly showed 
p-type conductivity in tunneling I(V) spectroscopy. The 
topographies were measured with a room temperature 
UHV-STM (P < 2 x 10 -11 torr) on an in situ cleaved 
(110) surface. 

Figure EJa) shows the electronic topography of one of 
the pairs in the ionized state (the other observed pairs 
are similar). In the ionized configuration, it is impossible 
to distinguish locations of the two dopants. The poten- 
tial from the double charge of the two ionized dopants 
induces an apparent round elevation 1.7 times larger 
than that of a single ionized Mn under the same imag- 
ing conditions pi . In the neutral configuration, how- 
ever the presence of two dopants can be clearly identified. 
Further two examples of close, clearly identifiable Mn-Mn 
pairs illustrate the resilience of the Mn wave-function to 
interaction with nearby Mn dopants [Fig. |2fc,d)]. 

A schematic model of the GaAs (110) surface is shown 
in Fig. Efb). The surface locations of the two Mn of 
the first pair [shown in Fig. [2f a,c,e,g,i) and separated by 
1.4 nm] are indicated by Mno and Mni; the Mn atoms 
themselves are located in the fifth sub-surface layer, and 
are well separated from neighboring Mn dopants. The 
other pair [shown in Fig. (2{d,f,h,j) and separated by 
0.8 nm] is indicated in figure Efb) by Mno and Mn2, and 
is likewise in the fifth sub-surface layer as well. This 
combination has the smallest separation of those Mn-Mn 
pairs that we were able to identify. These measurements 
show that in the neutral state, the wave- functions of the 
two Mn acceptors retain their crosslike shape even when 
they are separated by a distance smaller than the wave- 
function's effective Bohr radius ao « 0.9 nm. 

Figure Efe,f ) shows the topography calculated with the 
tight-binding model (TBM) |8|, |9( for the two pairs shown 
in Fig. Etc,d), and for Mn spins parallel to each other. 
The calculation is averaged over the orientation of the 
two parallel Mn spins relative to the crystal's lattice. 
Qualitatively the crosslike shape is clearly evident, and 
the agreement between the calculations and the measure- 
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merits is as good as found for a single Mn dopant in Ref.0. 
The TBM is based on the deep level model of Vogl and 
Baranowski[25j and is applied to a bulk-like Mn accep- 
tor. The dangling sp 3 -bonds from the nearest-neighbor 
As hybridize with the Mn <i-states of Tis character. The 
antibonding combination of these becomes the Mn ac- 
ceptor state. Coupling to the (i-states of Ti2 character is 
weak, and hence neglected. The hybridization strength 
is fully determined by the acceptor level energy. 

Comparison with theoretical calculations based on the 
TBM permits a quantitative evaluation of the effect of 
the Mn-Mn interaction on the measured topography. 
The topography for a "non-interacting baseline" is con- 
structed by adding together the topography of two single, 
isolated Mn displaced by the pair separation. This base- 
line is then subtracted from calculations of Mn pairs with 
differing spin orientations. Figure Efg,h) shows the differ- 
ence between Fig.|2{e,f) and the non-interacting baseline. 
The quantity A is the ratio of the largest difference shown 
in Fig. EJg,h) to the largest value of the topography for 
the non-interacting baseline, and even for the close pair 
is less than 1/3. In previous work the spectral and spa- 
tial differences between dopant spin pairs with parallel 
and antiparallel spins were predicted 0, [2^. Here the 
spectra could not be measured with sufficient resolution 
to distinguish the pair spin orientations, and the differ- 
ences expected between parallel and antiparallel spins are 
small. Figure EJij) shows the same as Fig.[2{g,h), except 
that the Mn spins are antiparallel. The differences be- 
tween parallel and antiparallel are of the order of 10%, 
which is not resolvable in our measurements. 

Now that the robustness of the anisotropic crosslike 
shape of the Mn hole wave-function has been clearly es- 
tablished, we explore the implications for spin-spin cou- 
pling mediated by these hole wave- functions in an impu- 
rity band. The experimental data acquired with STM is 
a two-dimensional slice along a (110) plane of the entire 
three-dimensional wave function. As a result, any estima- 
tion of the directional dependance of the wave-function 
overlap taken directly from the STM experiment would 
be incorrect. Instead we quantify the directionally de- 
pendent overlap of the wave-functions by calculating the 
bulk-like Mn-acceptor wave-function within a four-band 
Luttinger-Kohn envelope-function model (EFM) as well 
as the tight-binding model. The EFM uses the zero-range 
potential modeling, including a cubic correction as sug- 
gested in Ref. The ground state of the Mn acceptor 
can be approximated as four-fold degenerate with a to- 
tal momentum of the valence hole F = 3/2 and has the 
symmetry of the top of the valence band Tg j^l • We ne- 
glect possible effects caused by the presence of the (110) 
surface and quantum spin effects from the exchange in- 
teraction between the Mn 2+ 3d 5 core and the hole. 

The calculated radial dependance of the overlap of non- 
interacting Mn wave-functions for three crystallographic 
directions is presented in Fig. [31 The graph shows a 
nearly exponential decay of the overlap integral with sep- 
aration, but characterized by a directionally-dependent 
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FIG. 3: (Color online) Separation-dependent overlap of non- 
interacting Mn acceptors calculated for separations along 
three crystallographic directions using the envelope-function 
model (EFM) and the tight-binding model (TBM). 
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FIG. 4: (Color online) Directional dependance of the Mn 
wave-function overlap for Mn-Mn separations D calculated in 
the EFM. The curves are normalized to the maximum value 
N given in Table [0 



decay constant. Thus the anisotropy of the overlaps in- 
creases at larger separations between the Mn dopants. 
The calculated directional dependence of the overlaps of 
non-interacting Mn wave-functions for various Mn-Mn 
separations is shown in Fig. QJ The maximum of the 
overlap occurs when the impurities are located along the 
[111] direction, whereas the minimum occurs along the 
[001] direction. The TBM and EFM show similar qualita- 
tive behavior, however the results differ slightly in value. 
The EFM zero-range potential model underestimates the 
magnitude of the wave-function anisotropy compared to 
that observed in experiment and obtained with TBM. 

The wave- function overlap Qi( xyz ) for Mn-Mn pairs 
grown on (xyz) oriented substrates (Table [I]) is estimated 
by averaging the curves Ei for directions perpendicular to 
(xyz). These calculations suggest that the wave- function 
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TABLE I: Calculated values of the Mn wave- function overlap 
from the EFM. S is the curve number in figure 0] D is the 
Mn-Mn separation, N is the normalization coefficient, and 
^ xyz) is the averaged overlap integral for Mn pairs grown on 
a (xyz) substrate. 
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overlap on average is different for Mn (5-doped layers 
grown on differently oriented substrates. The density 
where the Mott metal-insulator transition occurs in an 
impurity band is determined by the overlap of localized 
wave functions. The anisotropic overlap of the Mn wave- 
functions will produce a directionally-dependent density 
threshold for percolation as well. Thus the critical con- 
centration for the metal-insulator transition will be lower 
for Mn J-doped layers grown on (111) or (110) substrates 
compared to ^-doped layers grown on (001) substrates. 
As the Curie temperatures of metallic GaMnAs are much 
higher than those of insulating GaMnAs, the Curie tem- 
peratures and other magnetic properties for J-doped lay- 



ers should be strongly dependent on the substrate orien- 
tation, with (111) substrates yielding higher Curie tem- 
peratures than (110) substrates or the currently-used 
(001) substrates. 

In conclusion, we have experimentally demonstrated 
that the crosslike shape of the Mn persists in groups of 
Mn with short Mn-Mn separation. This strongly sup- 
ports the picture of impurity-band conduction and spin- 
spin coupling at Mn doping densities corresponding to 
ferromagnetic GaMnAs. We suggest that the anisotropy 
of the Mn wave-function will substantially influence the 
carrier density of the Mott metal-insulator transition in 
Mn J-doped layers grown on differently oriented sub- 
strates. We expect that Mott transition will occur at 
lower Mn concentrations in layers grown on (111) sub- 
strates and at higher concentrations in layers grown on 
(001) substrates, leading to higher Curie temperatures 
for (lll)-grown than (OOl)-grown 5-doped layers. These 
results have broad implications for all acceptor-acceptor 
interactions in zincblende semiconductors, and especially 
for hole- mediated ferromagnetic semiconductors. 
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